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We examine the magnetic phase diagram of iron pnictides using a five band model. For the 
intermediate values of the interaction expected to hold in the iron pnictides, we find a metallic 
low moment state characterized by antiparallel orbital magnetic moments. The anisotropy of the 
interorbital hopping amplitudes is the key to understanding this low moment state. This state 
accounts for the small magnetization measured in undoped iron pnictides and leads to the strong 
exchange anisotropy found in neutron experiments. Orbital ordering is concomitant with magnetism 
and produces the large zx orbital weight seen at F in photoemission experiments. 
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One of the common features in most iron pnictides 
is the appearance of unusual antiferromagnetism (AF) 
in the undoped compounds. In this state the magnetic 
moments order with momentum Q = (7r,0), namely an- 
tiferromagnetically in the a;-direction and ferromagneti- 
cally in the y-direction[lj. The Curie temperatures are 
high Tjy 130 — 200K, while the measured magnetic 
moment m is small m ^ 0.3 — l.OfiB- A structural tran- 
sition at Ts > Tn accompanies magnetism and the sys- 
tem shows metallic behavior in the magnetic state. In 
spite of the current hot debate, the weak0, or strong 
coupling[3, ^ nature of magnetism is not clear yet. Ab- 
initio calculations [2] generally report magnetic moments 
™ > much larger than experimentally measured. 

Proposals to explain the small magnetic moment include 
frustration[5], weak order [61], antiphase boundaries jTi], op- 
posite orbital magnetizations [1, 1^ or the interplay be- 
tween frustrated and non- frustrated bands [loj. 

Unexpectedly, very anisotropic nearest neighbor ex- 
change constants Jy <IC J^, with Jy even slightly ferro- 
magnetic, have been necessary to describe neutron scat- 
tering results [l_lj]. Orbital ordering (00) was proposed 
early on within a Kugel-Khomskii description ll2j and ar- 
gued to be behind both the strong anisotropy |13 , [l4 and 
the structural transition 14, 1^. Recent experiments [lil 



17| have also been interpreted as manifestations of orbital 
ordering. In particular ARPES[l6l| experiments show 
that the Fermi surface at F has a predominant zx or- 
bital component. 

Here we study the magnetic phase diagram of a five 
band Hamiltonian at the mean field level. Different mag- 
netic regimes, mostly coexistent with 00, are found. At 
intermediate Hubbard interaction, the system is AF and 
metallic and shows two different magnetic phases. In par- 
ticular, we find a low moment (LM) phase which accounts 
for the small magnetization measured on iron pnictides. 
The LM arises as a consequence of partial cancellation of 
antiparallel orbital magnetic moments. The anisotropy of 
interorbital hoppings is the key to explain this phase. To 
compare with neutron results, we estimate the anisotropy 



of the exchange interactions finding it to be strong only 
in the low moment state. We also show that the 00 has 
a small contribution to the anisotropy of this state. On 
the other hand, the 00 enhances the zx orbital over the 
yz orbital weight in the Fermi surface. 

We start from the 5-band interacting Hamiltonian 
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Here i,i label the Fe sites in the Fe unit cell. 7 and /3 
refer to the five Fe-d orbitals yz, zx, xy, 3z^ — and 
x^ — included in the model, and a to the spin, x and y 
axis are directed along the Fe bonds. The kinetic energy 
term includes hopping up to second neighbors with the 
hopping arnplitudes calculated within the Slater-Koster 
frameworkjisj as detailed in Ref . . Both direct Fe-Fe 
and indirect (via As) hoppings determine the magnitude 
of the hopping amplitudes. This tight binding model[l3| 
gives good account of the band structure found in den- 
sity functional theory, including the orbital content of 
the bands, with a reduced number of fitting parameters. 
The hopping amplitudes depend on the angle a formed 
by the Fe-As bonds and the Fe-plane. In the follow- 
ing, a = 35.3° corresponding to the regular tetrahedra is 
used, except otherwise indicated. The interacting part of 
the Hamiltonian includes the intraorbital U and the in- 
terorbital U' interactions, as well as the Hund's coupling 
J and pair hopping J' terms. The pair hopping inter- 
action is written for completeness, but it does not enter 
in the mean-field approximation used below. Energies 
are given in units of {pda)'^/\ed — Cpl ^ 1 eV, with pda 
the a overlap between the Fe-d and As-p orbitals and 
|e<i — fpl their energy difference [l^. We assume that AF 
takes place with Q — (tt, 0) momentum, as experimen- 
tally observed, and treat the Hamiltonian at the mean 




FIG. 1: (Color online) Calculated (tt, 0)-mean field magnetic 
phase diagram of the five band Hamiltonian in Eq. [1] Pat- 
terned areas (for U > 3) correspond to gapped states. The 
following color code applies: gray corresponds to the para- 
magnetic (PM) state; red to a low magnetic moment (LM) 
state showing orbital magnetizations with opposite signs; 
blue to high magnetization (HM) with parallel m-y which cor- 
responds at large [/ to an Sz = 2 state; white is an interme- 
diate = 1 moment state; and green is the Sz ~ state. 
The Sz = 2, Sz ~ 1 and Sz ~ states are illustrated from 
top to bottom on the right. Magnetic moment m and gap are 
assumed to be finite when larger than 0.001. 



field level keeping only the spin and orbital-diagonal av- 
erage terms [20] 
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where k runs over the Fe Brillouin zone and cr = ±1. 
This corresponds to a magnetic moment m = "^^rriy, 
in units of /is. We assume that the relation U' = U — 



2J from rotational invariance[21[ holds and study the 
phase diagram in the J/U versus U space. Ferro-orbital 
ordering, denoted simply by 00 in the following, and 
AF with momentum Q are the only symmetry breaking 
states allowed in our calculation. We focus on undoped 
systems with n — — 6. 

Fig. [T] displays the phase diagram with the param- 
agnetic (PM) state shown in gray. Different magnetic 
phases with momentum Q are represented in green, red, 
blue and white and correspond to states with vanishing 
(5^ = 0), low (LM), high (HM) and intermediate, de- 
noted Sz — ^, magnetic moment, respectively. Shaded 
areas (for U > 3) correspond to insulating states with 
a gap at the Fermi level. The PM state, which shows 
no 00, survives up to a J-dependent critical value of U, 
Uc{J) ~ 1.8. The high value of Uc{J) is not expected in 
a nesting scenario and suggests a strong coupling origin 
of AF. Moreover, the magnetic states found for interme- 
diate values oi U ~ 2 — 3 are metallic. 

For large U the system evolves towards the atomic limit 
(see the sketches in Fig. [1]) with a well defined filling of 
orbital and spin states. For J ^ O.Olt/, an intermedi- 



FIG. 2: (Color online) (a) Total magnetic moment for three 
different values of J as a function of U. (b) Orbital magnetic 



moments m-y, (c) 00 



and gap, and (d) orbital fill- 



ing as a function of U for J/U = 0.07. The low magnetic 
moment state arising for 2 < U < 2.8 is due to the par- 
tial cancellation of opposite magnetic moments on different 
orbitals. 



ate magnetic moment Sz = I and insulating state ap- 
pears. In this state there is no 00 (7135,2 , n^2_y2 ~ 
2,nyz,nzx ~ 1) and myz,rnzx ~ l/^B- When J increases, 
a spin from one of the doubly occupied states is promoted 
to the xy orbital and we recover the 5^ = 2 atomic limit 
with 713,2 _y 2 2, nz2-r^,nyz,nzx,nxy 1 and very weak 
00 [see Fig. IHb)]. In this state all the spins in the half 
occupied states are parallel. This tendency at large U 
can be appreciated in Fig. [2] The HM state survives 
when reducing U for sufficiently large values of J/U. At 
intermediate values of U the system is metallic and shows 
00. From Fig. [2] (a) it is apparent that, within this HM 
state, achieving the small magnetization m < Ifis re- 
ported experimentally is only possible by fine-tuning, as 
in the weak order state discussed in Ref. @. 

The Sz — I and Sz — 2 atomic states have been 
widely used as starting points in previous works, where 
the Hund's rule is assumed to hold. As shown in Fig. [T] 
for large U and J ^ 0, an Sz — state with a strong and 



positive 00 Uyz » Uzx appears. In this state Hund's 
rule is violated. The electrons in the xy and yz orbitals 
have antiparallel spins, see sketch at the right. Most im- 
portantly, at intermediate U and not too large J /U , we 
find an LM phase (with 777 < \^b) with negative 00 
7iyz < Tizx ill which the individual orbital magnetizations 
have opposite signs. In this state 777j,2, m^2_y2 < 0, while 
777^2;, irixy > (7713^2 ~ 0). This is shown in Figl2] (b) 
for J/U = 0.07 within the range 1.8 < U < 2.7 (between 
the vertical dotted lines). The low or almost vanishing 
total magnetic moment in these two states arises due 
to the partial cancellation of otherwise relatively strong 
magnetic moments in close to half-filling orbitals. 

An AF solution with opposite orbital magnetizations, 
stabilized by the formation of large multipoles of the 
spin magnetization, was previously found in an ab-initio 
calculation Isl for LaOFeAs. Violation of Hund's rule in 
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FIG. 3: (Color online) (a) Dependence on a of the first near- 
est neighbors hopping amplitudes mostly responsible for the 
stability of the low-magnetic moment state and the orbital 
ordering. For other hopping amplitudes see [lH. (b) Direct 
Fe-Fe and indirect (via As) contributions to ty^^y^ and tyz,yz 
showing a strong anisotropy. (c) Sketch of the low magnetic 
moment state on three neighboring sites of the Fe plane and 
the largest hoppings involved in its stabilization. 



the pnictides has also been discussed within the context 
of a two-orbital Heisenberg niodel[9] as a consequence 
of large interorbital exchange J^,/?. In the latter work, 
isotropic J^^^ were used and the columnar AF order was 
not found. As discussed below, in our model the sta- 
bility of the Sz — and LM states can be understood 
within a strong coupling point of view by considering 
the anisotropy of the interorbital exchange interactions 
J^^p (xt'^^ which override Hund's rule. 

In an AF state with momentum Q = (tt, 0), each 
changes sign between nearest neighbors along the x direc- 
tion, but not along the y direction. As shown in Fig.[31[a), 
the hopping between yz and x'^ — y^ along the x direction 
\iyz 2:2-2/2 1 is large, while it vanishes along the y direction. 
Therefore, parallel magnetic order of yz and — y'^ is 
favored as the system gains AF energy along x without 
any cost in energy associated with the ferromagnetic or- 
dering along y. On the contrary, the hopping between 
yz and xy is large along the y direction {f^z^xyi van- 
ishes along X. This implies that parallel ordering of the 
yz and xy magnetic moments would cost exchange en- 
ergy along y, while there would be no gain whatsoever 
along X. A configuration with opposite signs of myz and 
ruxy shows relative AF ordering of these two magnetic 
moments along y (see sketch in Fig. [3]) saving exchange 
interaction. From symmetry, \t^ i_ 'i\ = \t^.. „2 .,2! and 



zx,xy I 



\^z,xy \ which favors antiparallel orientation be- 
tween rUzx and rUx^-y^ and parallel orientation between 
rUzx and rUxy Interactions involving — are more 
frustrated, n^z'^_j.2 ^ 2 and therefore its magnetization 
is small in the LM state. At large J/t/, Hund's rule 
dominates and the lowest energy correspond to the HM 
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FIG. 4: (Color onUne) (a) Anisotropy J^ff/J^s of the mag- 
netic exchange from a strong coupling approach and (b) or- 
bital ordering function of U and J. The 
anisotropy appears mostly within the Sz = and low mo- 
ment states (green and red (LM) regions in Fig[T|). The or- 
bital ordering accompanies the magnetization within a wide 
range of parameters but is not correlated with the exchange 
anisotropy. 



state. 

Negative 00 appears both in the LM and HM regimes 
[see Fig. IDJb)]. This produces a splitting of zx and yz 
bands at F leading to larger weight of zx at the Fermi sur- 
face (see supplementary material[22'|) in agreement with 
ARPES experiments[16]. 

As discussed previously by Lee, Yin and Ku [ll] based 
on ab-initio calculations, 00 originates in the anisotropy 
of the yz (zx) intraorbital first nearest neighbor hop- 
ping which favors a large (small) magnetic moment in the 
yz {zx) orbital. Unexpectedly, t^^ y^ > ty^ y^ {tl^^^zx > 
^zx zx)- This relationship is opposite to the one used in 
early proposals of 00 in iron pnictides 12,1^. Our tight- 
binding modell^ does not only reproduce Lee, Yin and 
Ku [IJ] results, but it also allows to understand its origin. 
From Fig.[3l[b) it can be seen that the small value of t^z^yz 
{— ^zx zx) comes from the cancellation of large direct Fe- 
Fe and indirect (via As) contributions with opposite sign, 
which add with the same sign for ty^ y^ {— t^^ .,^) 



In order to compare with neutron experiments ll| we 
make connection with a Heisenberg model and estimate 
the anisotropy of the exchange interactions. Within a 
strong coupling approach, J^'J oc (fi^'p^, and assuming 

JcS = J2-y,p Jj'}^^,i^i3,j/{SiSj) with s^j and Sj the or- 
bital and total spin at site i, we have calculated the ratio 
J^ff/Jgff. The results for the anisotropy are shown as a 
function of U and J/U in Fig. |4] (a). Large anisotropy 
J^ff ^ J^g, including negative values of J^g/J^^, is only 
found in the LM state. A value comparable to the 
experimental ll[ Jgg/Jgg ~ —0.11 is found in the LM 
state close to the transition to the HM state. In the 
HM state, the largest calculated anisotropy corresponds 
to J^jj/ J^g ~ 0.88, very far from the experimental value. 
Comparing with Fig. |4] (b), it is apparent that the ex- 
change anisotropy in Fig. 0] (a) is not strongly related to 
00 (which appears in most of the phase diagram), but 
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to the presence of antiparallel orbital magnetizations, ex- 
cept for the insulating strong coupling solution Sz = 
where both effects are present. This is easily understood 
as the contribution of products s-^s^ with positive and 
negative sign cancel when summed up. Following a sim- 
ilar procedure we have also estimated the ratio of the 
exchange between second and first neighbors J^^ / J^^. 
Its value, not shown, is around 1/2 in HM state but de- 
creases in the LM phase. This is also consistent with 
the experimental [m j'^ / J^g ^ 0.38. This suggests that 
second nearest neighbors are important to stabilize the 
HM, but not the LM state. 

Several works have emphasized the influence of the 
Fe-As angle a on determining the properties of iron 
pnictides[23|. We have calculated the phase diagram 
for a — 37.2° (elongated tetrahedra) and for a — 29.9° 
(squashed tetrahedra as found in LaFePO). For the elon- 
gated case the phase diagram is very similar to the one 
discussed here for the regular tetrahedra. On the other 
hand, for a = 29.9° the low magnetic moment is less 
stable and reduced to a very small portion of the phase 
diagram. It is interesting to note that AF is absent in 
LaFePO. The worse stability of the LM state for squashed 
tetrahedra can be understood by looking at the angle 
dependence of the hopping parameters (see Fig. Ela)): 
\'t%z,xy\ which helps stabihze the LM phase strongly de- 
creases with decreasing a. 

In conclusion, we have studied the mean field mag- 
netic phase diagram of a five orbital model for the iron 
pnictides. Several magnetic phases appear for different 
values of Hubbard and Hund interactions. 00 is present 
for a wide range of parameters but does not seem to have 
a strong effect on the anisotropy of the exchange inter- 
actions in the metallic region. A metallic low magneti- 
zation state with antiparallel orbital magnetic moments 
is found for intermediate values of U . It is stabilized by 
the anisotropy of the interorbital hoppings, ultimately re- 
lated to the symmetry of the orbitals and the tetrahedral 
coordination of the As atoms. This state is consistent 
with the measured small magnetic moment [J, the large 
weight of zx found in the Fermi Surface around r[16|, 
and the anisotropy of the exchange interactions [ll|. 

Our results uncover antiparallel orbital magnetizations 
as a new source of anisotropy connecting not previously 
related experiments: the anisotropy measured with neu- 
tron scattering and the observed low magnetic moment. 
They suggest a strongly correlated origin of the mag- 
netic state different from both the nesting scenario and 
the usual Heisenberg description in terms of the atomic 
moment and point to the need of describing correctly 
the individual exchange interactions between the orbital 
magnetic moments. They also stress the importance of 
including all five Fe d-orbitals in the description of the 
pnictides and question the validity of band models which 
mimic the band structure and Fermi surface with less or- 
bitals as the hopping amplitudes are expected to differ 



considerably from the correct ones. 

To the best of our knowledge, we are not aware of any 
technique that can make a straightforward direct mea- 
surement of the low magnetic moment state that we find. 
However, we expect it to show up indirectly as we have 
already seen in the anisotropy of the exchange interaction 
observed with neutron experiments. 
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